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Abstract 

Universal Bernoulli polynomials are introduced and their number- 
theoretical properties discussed. Generalized Euler polynomials are also 
proposed. 

1 Introduction 

The purpose of this paper is to introduce new classes of polynomial sequences 
which represent a natural generalization of the classical Bernoulli and Euler 
polynomials. In particular, we introduce the notion of universal Bernoulli poly- 
nomials. 

The Bernoulli polynomials [5] are defined by the generating function 



1 jfe! y ' 



k=0 

For x = 0, formula reduces to the generating function of the Bernoulli 

numbers: 

^ =£-!''■ <'- 2 > 

fe=0 

The Bernoulli polynomials are also determined by the two properties 

D B n (x) = nB n _\ (x) , (1.3) 

AB n (x) = nx n -\ (1.4) 

with the condition Bq (x) = 1. Here Dis the continuous derivative, A := (T — 1) 
is the discrete one, and T is the shift operator, defined by Tf (x) = f (x + 1). 
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Standard references are the monographs JS| an d 03 an d tne books [23 and 
|27|. By virtue of relation (|1.3|l . the Bernoulli polynomials belong to the class 
of Appell polynomials. Formula l|1.4|) is the basis for the application of these 
polynomials in interpolation theory |36|. They also play an important role in 
the theory of distributions in p-adic analysis (2H|, in the theory of modular 
forms in the study of polynomial expansions of analytic functions UJ, etc. 
Several generalizations of them have also been proposed. Particularly important 
is that of Leopold and Iwasawa, motivated by a connection with the theory of p— 
adic L-functions [2E]- Another example is provided by the work of Carlitz 
Recently, new applications of the Bernoulli polynomials have also been found in 
mathematical physics, in connection with the theory of the Korteweg-de Vries 
equation [23 and Lame equation |2"T . and in the study of vertex algebras [17]. 

The Bernoulli numbers are relevant in several branches of number theory, in 
particular to compute rational values of the Riemann zeta function |14j . |2()| . 
in the theory of cyclotomic fields [45| and, since Kummer's work, in connection 
with Fermat's last Theorem pTJ]- They are also useful in singularity theory |3] 
and in connection with Coxeter groups 0]. Standard applications in algebraic 
topology are found in the computation of Todd characteristic classes and in the 
Hirzebruch signature Theorem, as well as, more recently, in complex homology 
theory [S], [S7]. In the last years Bernoulli number identities have found applica- 
tions in Quantum Field Theory J§| and in the computation of Gromow-Witten 
invariants |23| . 

Inspired by the important works by Clarke JS], Ray [33 and Adelberg 
PQ on universal Bernoulli numbers, in this paper we introduce the following 
generalization of the Bernoulli polynomials. 

Definition 1. Let ci, C2, ... be indeterminates over Q and let 

F( S ) =S + Cl y +C 2 y + ... (1.5) 

Let G (t) be the compositional inverse series: 

G(t) = t~ Cl j + (3cl-2c 2 )j + ... (1.6) 

so that F (G (t)) = t. The universal higher-order Bernoulli polynomials 
B k,a ( x > c i> •••> c n) = B k,a ( x ) are defined by 

(4)>=x>?.(4- <"> 

v w/ fc>0 

where x, a G M. 

We immediately observe that, when a — 1, a — (— 1)*, then F (s) — 
log (1 + s) and G (t) = e — 1, and the universal Bernoulli polynomials and num- 
bers reduce to the standard ones. For a 6 Z, a = (—1)' we reobtain the 
higher-order Bernoulli polynomials, which have also been extensively studied 
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(e.g. in EH], [12, [13], [H]). When x = 0, a £ R and q = (-1)' formula (|T7jl 
reduces to the standard Norlund polynomials, and for a € Q to the universal 
Norlund polynomials. 

By construction, the polynomials B%i (0) = Bk € Q [ci, c„], where Bk 
are the universal Bernoulli numbers introduced by Clarke in |15) . The name 
comes from the fact that G (F (si) + F (S2)) is the universal formal group 
the series F is called the formal group logarithm and G the formal group expo- 
nential. The universal formal group is defined minimally over the Lazard ring 
L, which is the subring of Q [ci, c„] generated by the coefficients of the power 
series G (F (si) + F (52))- The universal Bernoulli numbers have been recently 
investigated in the context of algebraic topology, in particular in complex cobor- 
dism theory, in US] and [2Zj, where the coefficients c„ are identified with 
the cobordism classes of CP". They also obey generalizations of famous con- 
gruences valid for the classical Bernoulli numbers, like the celebrated Kummcr 
and Clausen- von Staudt congruences 0, |27) . 

We will show that the universal polynomials l|1.7fl possess several interesting 
properties. An intriguing one is a generalization of the Almkvist-Meurman [2] 
and Bartz-Rutkowski congruences, before known only in the case of classical 
Bernoulli polynomials. 

In this paper we will study in particular a class of polynomials of the type 
l|1.7fl constructed from the finite operator calculus, as formulated by G.C. Rota 
and S. Roman 0U|. For a recent review of the vast literature existing on 
this approach (also known, in its earlier formulations, as the Umbral Calculus), 
see e.g. ^H]- The Bernoulli-type polynomials of order q considered here corre- 
spond to a class of formal power series G (t) obtained from the difference delta 
operator of order q (denoted by A q ) introduced in [^j , pH] • These polynomi- 
als, parametrized by a real variable a, are uniquely determined by the relations 

DB* ta (x)=nB*_ lta (x), (1.8) 

A q Bl a (x)=nx n - 1 , (1.9) 

and by the condition Bq (x) = 1. Generating functions as well as many com- 
binatorial identities will be presented. For q = 1, a = 1, they coincide with the 
standard Bernoulli polynomials. In the same spirit, a generalization of Bernoulli 
polynomials of the second kind is also proposed. Although the use of Roman- 
Rota's formalism, strictly speaking, is not indispensable for formulating our 
results, it has the advantage of providing a natural and elegant language which 
allows the mathematical treatment of many polynomial sequences to be unified. 

The original motivation for studying the polynomials proposed here is dis- 
crete mathematics. Indeed, they satisfy certain linear difference equations of 
order q in one variable, defined in a two-dimensional space of parameters. In 
|53] . [33 a version of Rota's operator approach, based on the theory of rep- 
resentations of the Heisenberg-Weyl algebra introduced in (32], 021) nas been 
outlined. One of the goals of this paper is to establish a connection between 
Rota's approach and the theory of Appell and Shcffcr polynomials. The basic 
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sequences associated with the delta operators A q are studied in detail and some 
of their combinatorial properties derived. Their connection constants with the 
basic sequence {x n } (and the associated inverse relations) define a class of gener- 
alized Stirling numbers of the first and second kind. The problem of classifying 
all Sheffer sequences associated with these delta operators is essentially open. 

Euler polynomials and numbers (introduced by Euler in 1740) also possess 
an extensive literature and several interesting applications in Number Theory 
(see, for instance, [T^], ESj EH)- Prom many respects, they are closely related 
to the theory of Bernoulli polynomials and numbers. A generalization of Euler 
polynomials is also presented. 



Definition 2. The Euler-type polynomials are the Appell sequence generated 
by 



E 

fc=0 



(1.10) 



where a, lu G 



Definition 3. The Euler 



ibers E 



(0) are defined by 



k=0 



k\ 



(0)^ 



(1.11) 



For q = a = 1, to = 2, we reobtain the classical Euler polynomials. It will 
be shown that these new polynomials possess many of the properties of their 
classical analogues. 

The paper is organized as follows. In Section 2, some properties of universal 
Bernoulli polynomials are discussed. In Section 3, a brief introduction to finite 
operator calculus is presented, with a discussion in Section 4 of the delta opera- 
tor theory. In Section 5, generalized Stirling numbers are considered. In Section 
6, new Bernoulli-type polynomials of the first kind are introduced. They are 
used in Section 7 to generate sequences of integer numbers. The Bernoulli-type 
polynomials of the second kind are introduced and discussed in Section 8. The 
proposed generalization of Euler polynomials is discussed in the final Section 9. 



2 Some properties of universal Bernoulli poly- 
nomials and numbers 

The universal Bernoulli polynomials l|1.7f) possess many remarkable properties. 
By construction, for any choice of the sequence {c n } ngN they represent a class 
of Appell polynomials. 

We easily deduce the following results. The Appell property is expressed by 
the identity: 

B% a (x + V ) = J2 ( I ) B l* (*) y n ~ k - (2- 1 ) 



4 



In particular, for x = we obtain an useful characterization of B„ a (x). 

Lemma 1. The higher-order universal Bernoulli polynomials are expressed in 
terms of the higher-order universal Bernoulli numbers by the relation 

k=0 ^ ' 

From relation (|2.1|1 . applying the operator {jjij^j we get 

{x + y) = J2 ( I ) B l« (*) B n-k, b (V) ■ (2-3) 



fc=0 



and, for (a + b) = 0, 

(* + 2/)" = E ( I ) B la (*) (W) 

t — n V / 



(2.4) 



A recurrence formula is valid 



G(t) 

The following multiplication property holds: 

BV{ax) = -£^BV{x). (2.6) 

In the case of classical Bernoulli polynomials it reduces to the famous Raabe 
multiplication Theorem. 

The values B% \ (0) := B n correspond to the well-known universal Bernoulli 
numbers. Here we mention only two of their most relevant properties. 

i) The universal Von Staudt's congruence [TK] . 

If n is even, 

c «/(p-i) 

B n =- Y mod Z[ Cl ,c 2 ,...]; (2.7) 

p 

p-l|n 
pprimc 

If n is odd and greater than 1, 



„n I Jl—3- 

B n = 1 ^ 3 mod Z[ Cl ,c 2 ,...]. (2.8) 



When c„ = (—1)", the celebrated Clausen- Von Staudt congruence for Bernoulli 
numbers is obtained. 

ii) The universal Kummer congruences p^. The numerators of the classical 
Bernoulli numbers play a special role, due to the Kummer congruences and to 
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the notion of regular prime numbers, introduced in connection with the Last 
Fermat Theorem (see, for instance, |27|). The relevance of Rummer's congru- 
ences in algebraic geometry has been enlightened in |S] . More general versions of 
these congruences for the classical Bernoulli numbers are known in the literature 
|44| . As shown by Adelberg, the numbers B n satisfy an universal congruence. 
Suppose that n^0,l (mod p — 1). Then 

'" ' - ™-c p _i mod pZ p [ci,C2,...] . (2.9) 



n + p — 1 n 

Almkvist and Meurman in j2j discovered a remarkable congruence for the 
Bernoulli polynomials: 

k n B n (~) e Z, (2.10) 

where k,h,n S N. Bartz and Rutkowski have proved a Theorem which com- 
bines the results of Almkvist and Meurman and Clausen- von Staudt. The 
main result of this section is a simple generalization of the Almkvist -Meurman 
and Bartz-Rutkowski Theorems for a class of Bernoulli-type polynomials. 

Theorem 1. Let h > 0, k > 0, n be integers and p > 2 be a prime number. 
Consider the polynomials defined by 



k 



v 1 k>0 

where G (t) is given by H.6]) . with c v -\ = 1 mod p and the other ci € Q. Then 

k n Bv(j) eZ[ci,C2,...], (2.11) 



where ( x ) = B v n (x) - B n . 

The following Lemmas are useful in the proof of Theorem ^ 
Lemma 2. If s is a positive integer, then 

E S- ( r 1, /P^[ci,C2,..]. (2.12) 

p\s 

It is a consequence of Fermat 's little Theorem. 
Lemma 3. Let n be a positive integer and p > 2 a prime number. Denote 

K 

x ^ / n 



fe=i 



p-1 



where [a] is the unique integer such that [a] < a < [a] + 1 . Then S p (n) = 
mod p if p — 1 \n and S p (n) = 1 mod p if p — 1 | n. 
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This Lemma is a particular case of an old result of Jenkins (see v °l- I; 
pag271). 

Proof of Theorem ^ It is a straightforward generalization of the proof 
of the Almkvist-Meurman Theorem given in [7] . The argument is based on the 
induction principle on t with s fixed. First observe that, by virtue of (|2.1|) for 

x — j and y — ^ we get 



t + 1 



E 

rn — 



m 



s m BZ 



For t = , the Theorem is obviously true. Assume that it holds for t. We 
obtain: 



s n BU 



E 

m— 



t + i 



n 
rn 



E 

m=0 



n 
rn 



s n B n 



n-l 

E 

m=0 



n 
rn 



6 JD Tt 



(2.13) 



Due to the induction hypothesis, we have only to prove that the congruence 
<|2.11[l holds for the second summation in (|2.13|) . 

Notice that under our hypotheses, for every odd m > 1 we have 



n 
m 



s m B m eZ[c 1 ,c 2 ,...}. 



•\n— l 

■/m>lodd 

Consider first n even. We shall distinguish two subcases. 

Let s be an even integer. Taking into account Lemma[21 we get 



n-l 

E 



m = 

m— 1 or m even 



n 
rn 



E^ M) /p=E E 

p—l\m 
^ n — 1 

— H + 



m/(p-l) 



p-i|, 



m=0 

p— l|m 



n 
m 



e z, 



where H € Z. Here Lemma 01 has been used for A' = ^— j- . Recalling that 

congruences ()2.7JI hold mod Z [ci, C2, ...], the Theorem is true for t + 1. A very 
similar reasoning applies when s is odd, considering separately the cases p > 2 
and p = 2. Now 



E 



m = 

m— 1 or m even 



n 
rn 



n-l 



ECr i3 /P=Ej E 



m/(p-l) 



,-l|rr 



,>\s 



m = 
i = l or m c\ 

p— l|m 



n 
rn 



n-l 



P f 3 m=0 
p>2 p— l|m 



m 



n + 1 



-ci 



"- 1 „m 

St 

m— 



n 
rn 



€ Z. 
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The use of mathematical induction completes the proof. The case n odd is 
very similar and it is left to the reader. The next result is a direct generalization 
of the Bartz-Rutkowski Theorem. 



Corollary 1. Under the hypotheses of Theorem^ we have that: 
if n is even 



(h\ n/p—l 




Strictly related to the AM Theorem is an interesting problem proposed by 
A. Granville. 

Problem: To classify the space of the polynomial sequences satisfying the 
Almkvist-Meurman property 



and to find a convenient basis. 

Although we will not address this problem explicitly we observe that The- 
orem Q provides a tool for generating sequences of polynomials satisfying the 
property (|2.14|) . More generally in Section 6 polynomial sequences satisfying 
Theorem ^ will be explicitly constructed using the finite operator theory. 

Many other congruences concerning the universal Bernoulli polynomials and 
numbers can be constructed. They will be discussed in a future publication. 

3 Finite operator calculus 

In this Section some basic results of the theory of finite difference operators 
and their relation with classical polynomial sequences are reviewed. For further 
details and proofs, see the monographs 0U|, (30j, where an extensive and modern 
treatment of this topic is proposed. 

Let T denote the algebra of formal power series in one variable t, endowed 
with the operations of sum and multiplication of series. Let V be the algebra of 
polynomials in one variable x and V* the vector space of all linear functionals 
on V . If L £ V* i following Dirac we will denote the action of L on p (x) £ V by 
(L\ p (x)). A remarkable fact of the finite operator calculus is that any element 
of .Fcan play a threefold role: It can be regarded as a formal power series, as a 
linear functional on V and also as a linear operator oriV. To prove this, let us 
first notice that, given a formal power series 




(2-14) 



OO 




(3.1) 
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we can associate it with a linear functional via the correspondence 

(/ (t)\ x n ) = a n . (3.2) 



Therefore 



In particular, since 



m _£<m*lt. ( 3. 3) 



fc=0 



(t k \ x n ) =n\S nik , (3.4) 
it follows that, for any polynomial p (x) 

(t k \ P (x))=p k (0). (3.5) 

where p k (0) denotes the fc-derivative of p(x) evaluated at x = 0. It is easily 
shown that any linear functional L € V* is of the form (|3.3I) . 

If we interpret t k as the k-th order derivative operator on V, given a poly- 
nomial p (x) we have that t k p (x) — p^ (x) . Hence the formal series 

oo k 
fe=0 

is also regarded as a linear operator acting on V . Depending on the context, t 
will play the role of a formal variable or that of a derivative operator. 

Now, some basic Definitions and Theorems of finite operator theory, neces- 
sary in the subsequent considerations, are in order. 

Definition 4. An operator S commuting with the shift-operator 

[S, T] = 

is said to be shift-invariant. 

Relevant examples of operators belonging to this class are provided by the 
delta operators. 

Definition 5. A delta operator Q is a shift-invariant operator such that Qx = 
const 0. 

As has been proved in |40| . there is an isomorphism between the ring of 
formal power series in a variable t and the ring of shift-invariant operators, 
carrying 

/W = Ey^Ef. (3.7) 

k=0 ' k=0 

In the following let us denote by p n (x) a polynomial of order n in x. 
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Definition 6. A sequence of basic polynomials for a delta operator Q is a 
polynomial sequence satisfying the following conditions: 



For each delta operator there exists a unique sequence of associated basic 
polynomials. 

Definition 7. A polynomial sequence is called a set of Sheffer polynomials for 
the delta operator Q if 

1) s (x) = c ^ 0; 

2) Qs n (x) = ns n -i (x). 

Definition 8. An Appell sequence of polynomials is a Sheffer set for the delta 
operator D . 

Any shift invariant operator S can be expanded into a formal power series 
in terms of a delta operator Q: 



with ak = [Spk (x)] \ x =q, where pk is the basic polynomial of order k associated 
with Q. By using the isomorphism l|3.7fl . a formal power series s (t) is defined, 
which is called the indicator of S. 

Remark 1. A shift invariant operator is a delta operator if and only if it 
corresponds, under the isomorphism i)3. 7|) . to a formal power series G (t) such 
that G (0) = and G 1 (0) ^ 0. This series admits a unique compositional 
inverse. 

The umbral formalism also allows us to characterize the generating functions 
of classical polynomial sequences of Sheffer and Appell type. 

If s n (x) is a Sheffer set for the operator Q, then there exists an invertible 
shift-invariant operator S such that 



where p n (x) is a basic set for Q. 

Let us denote by s (t) and q (t) the indicators of the operators S and Q. We 
will say that s n (x) is the Sheffer sequence associated with (s (t) , q (t)). The 
following result holds: 



1) po (x) = 1; 

2) p n (0) = 0V n> 0; 

3) Qp n (x) = npn-i (x 




fe>0 



S S n (x) = p n (x) , 




(3.8) 



where q 1 (t) denotes the compositional inverse of q(t). 
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Another characterization of a Sheffer sequence is provided by the identity 

Sn (X + y) = f \ ) S k 0) Qn-k (j/) ■ 

fc>0 ^ ' 

In the specific case of the Appell sequences, which obey the equation 

D s n (x) = nsn-i (x) , (3.9) 

there exists an invertible operator g (t) such that 

g (t) s n (x) = x n . (3.10) 

Therefore, we will say that s n (x) is the Appell sequence associated with g (t). 
Its generating function is defined by 

(3.11) 



The Appell identity is 



Sn ( x + y) = J2( l )sk{y)x n - k . (3.12) 



k=0 

We also recall the polynomial expansion Theorem. If s n (x) is a Sheffer set 
for the operators (s (t) ,q(t)), then for any polynomial p (x) we have: 

( S (t)q(t) k \ P (x)) 

= £ ij (*) ■ ( 3 - 13 ) 

fe>0 

In particular, if q n (x) is a basic sequences, this formula reduces to 

P( X )=Y,~ y -»(*)■ (3.14) 



k>0 



4 Difference delta operators 

In this paper, we will consider a specific class of difference operators, introduced 
in |33j . having the general form 

Q = A q = - ^2 a kT k , l,meZ, I < m, m-l = q, (4.1) 



a k=l 



where and a are constants. In order to fulfil the Definition of delta operator, 
we must assume that 

m m 

a k = 0, kak = c ' ( 4 - 2 ) 

k=l k=l 



11 



We also require that in the continuous limit A reproduces the standard deriva- 
tive -D; this implies c = 1, i.e. 



J2ka k = l. (4.3) 



k=l 

Definition 9. A delta operator of order q = m — I is a difference operator of 
the form \4-l\j which satisfies eqs. an d 14 -ty - 

We observe that eq. I|4.1|l involves m — I + 1 constants <Zfc, subject to two 
conditions l|4.2[l and (|4.3|) . To fix all constants a k we have to impose m — I — 
1 further conditions. A possible choice is, for instance 

m 

j q = J2k q a k = 0, q = 2,3,...,m-l. (4.4) 

fc=Z 

When conditions (|4.2I) . (|4.3|l and (|4.4|) are satisfied, the operator 1)4. 1|) provides 
an approximation of order p of the continuous derivative D, since 

rn—l m 
v ; k=l 

From now on, we will put a = 1. The Pincherle derivative of a delta operator 
is defined by the relation 

Q' = [Q,x}. 

Now, let us introduce a shift-invariant operator [3 such that 

[A,x0\ = l. (4.5) 

It follows that (3 = (A') -1 |H2] _ !!II]- Such an operator is invertible 0O| and 
finite. Indeed, if A is of order n, using the identity = 1, the action of (3 

on a monomial of order n is easily seen to be 

(3x n = x n - a^x"- 1 - 3 , (4.6) 
3=0 

where a™ are defined via the recursion relation 

i-i 

n 
.7 + 1 



1=0 



with 7 - = Y^k=i ^ a k- We see that (3 preserves polynomial structures. When 
A = D, j3 = 1. Other specific cases are listed below: for 

A = A+ = T-1, f3 = T- 1 ; (4.7) 
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for 

A = A~ = 1-T~ 1 , j3=T. (4.8) 

When 

T _ T -i fT + T~" 



A 2 = A* = ^— , /3=^-2— ) • («) 
Other nontrivial examples of higher-order operators, for instance, are provided 

by 

A 3 = - (T 2 - 2T + T- 1 ) , (4.10) 

A 4 = T 2 - ^T- 1 -T~ 2 , (4.11) 

A 5 =T 3 -2T 2 + 2T -2T- 1 +T~ 2 , (4.12) 

A 6 =T 3 ~2T 2 + 2T -T- 1 -T~ 2 + T~ 3 : (4.13) 

A 7 = T 4 - T 3 + T 2 - 2T + T- 1 - T~ 2 + T~ 3 . (4.14) 

Finally, we recall that Q is a delta operator if and only if it is of the form 
Q = DP, where P is an invertible shift-invariant operator It follows that 
the sequence of basic polynomials for Q is expressible in the form 

Pn (x) = xP- n x n - 1 . (4.15) 

5 A generalization of the Stirling numbers 

According to Rota's approach, with any delta operator of the form 1)4.10 — l|4.3fl it 
is possible to associate a sequence of basic polynomials. In perfect analogy with 
the classical theory of Stirling numbers, this allows to study generalizations of 
Stirling numbers. In this section, several concrete examples of these basic poly- 
nomials as well as generalized Stirling numbers are obtained. This construction 
represents a natural realization of the very general scheme proposed by Ray 
in |37| . who first studied the basic sequences and defined the universal Stirling 
numbers for any delta operator constructed via the complex homology theory. 
Here we will follow a different but equivalent philosophy: the basic sequences 
for the operators H4.1f) - H4.3[) are simply derived from the " discrete" representa- 
tions of the Heisenberg-Weyl algebra; the associated Stirling numbers furnish 
the connection coefficients between the basic sequences and the standard power 
sequence {x™}„ eN - 

We remind that the classical Stirling numbers of the first and second kind 
s (n, k) and S (n, k) are defined by the relations |38j 



(x) n = Y i 8(n,k)x* (5.1) 



13 



and 

n 

x n = J2S(n,k)(x) k , (5.2) 



fe=0 



respectively. Here (x) = x (x — 1) ... (x — n + 1) denotes the lower factorial 
polynomial of order n. The numbers S (n, k) also admit a representation in 
terms of a generating function: 



x n (e x - l) k 



As is well known, the lower factorial polynomials are the basic sequence associ- 
ated with the forward derivative: A {x) = n (x) -,. More generally, the basic 
sequence associated with the operator A q , denoted by (x)£, is [33| 

{x)% = {xP) n l. (5.3) 

Therefore, 

A(x)l=n(x)l_ 1 . (5.4) 
Explicitly, the first basic polynomials are expressed by 

(a;)g = (asjS) 1 = 1, 

(x)\ = (x/3) 1 1 = X, 
[x)l = (x(3) 2 1 = x 2 - 7 2 x, 
(a;)§ = {xpf 1 = x 3 - 3 l2 x 2 - ( 73 - 3 7 |) x, 



(x)l = (x(3) 4 l = x A - 6 72 .t 3 + (-4 73 + 15 7 |) x 2 + (-74 + IO7273 ~ x, 

(5.5) 

and so on, with 7 - defined by 

m 

7j =5^o fc ^, 7o = 0, 7l =l, j — 0, 1, 2, .... (5.6) 

fc=; 

Definition 10. TTie generalized Stirling numbers of the first kind and order q 
associated with the operators denoted by s q (n,k), are defined by 

a 

( Xn y = J2 s "(n,k)x k . ( 5J ) 

fc=0 

Remark 2. In the subsequent considerations, with an abuse of notation, we 
will use Roman numbers for denoting the values of q, when q appears as an 
index in sequences of polynomials or numbers. This to avoid possible confusion 
with exponents. 
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Observe that {x) n — (x) n and s 1 (n, k) — s (n, k). Since for any polynomial 
p (x) the following expansion holds: 



fe^O 

it emerges that 



s 



>,fc) = -U fc |«>. (5.9) 



A:! 



Definition 11. The generalized Stirling numbers of the second kind and order 
q, denoted by S q (n, k), are defined by 

n 

x n =J2S q (n,k)(x)l (5.10) 

k=0 

They admit the generating function 

E**(n,^ = ^f£, Ml) 

n—k 

where A q (t) = ■ aje J * is the indicator of the operator A g . As a consequence 
of formula (|3.14(l . one immediately obtains that 

S^n,k) = ^\^a je A \x n ^, m-l = q. (5.12) 
The polynomials (x n ) q satisfy the binomial identity 

+ = I )wiiy)l-k- (5-i3) 

k=0 ^ ' 

A relevant feature of relations l|5.7|l and (|5.1U|I is that they are inverse to each 
other. This immediately implies that 

8 mn = s q (n, k) S q (k, m) = Y,S q (n, k) s q (k, m) . (5.14) 

k k 

Remark 3. Due to formulae (|4.5|l and (|5.4|) . it is natural to interpret the 
operators A and x(3 as annihilation and creation operators acting on a finitely 
generated space, as noticed in 021 > |43j . For instance, 

(z/3)"-l=(at) n |0) = |n>. 
Ag(xf3) n -l = n\n-l). 
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6 Bernoulli— type polynomials and numbers of 
the first kind 



We propose the following Definition. 

Definition 12. The higher order Bernoulli-type polynomials of the first kind 
are the polynomials defined by the relation 

°° B q (x) 

£-*£r^ = (6.1) 

k=0 

with 

where Ylj—i o-j^ *s the indicator of a delta operator of order q. 

Remark 4. We will restrict to the case aj S Q, in order to ensure the validity 
of the results stated in Section 2. However, more general operator structures 
with real coefficients aj can be considered within the same philosophy. Observe 
that, when q = a = 1 we reobtain the classical Bernoulli polynomials. In the 
case a ^ 1 , q = 1 the polynomials l|6.2[l reduce to the classical higher-order 
Bernoulli polynomials. 

Definition 13. The Bernoulli-type numbers B\ are defined, for every q and 
a, by the relation 

oo B q 

E^p = J ?w a - ( 6 - 3 ) 

fc=0 

Remark 5. The operators J q (t) and the corresponding numbers for a = 1 
have been studied in |37] in the context of homological theories. Here again the 
construction of these operators is realized using a difference operator approach. 

We can now prove a result stated in the Introduction. 

Lemma 4. The polynomials (x) are uniquely determined by the two prop- 
erties 

DB" m (x) = mBl^ (x) , (6.4) 



A q Bl(x)=mx m -\ (6.5) 
with the condition Bq (x) = 1. 

Proof. We have 

B?n(x + h) = J2 7« (*) = £ ( v ) ^ 

!/=0 ' v=0 ^ ' 
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= E( y )h m -»Bl(x), (6.6) 

after having replaced v with m — v. Eq. (|6.5|l involves Bernoulli-type polyno- 
mials evaluated at shifted points. Substituting eq. (|6.6|l into eq. (|6.5|l and using 
the initial condition one can get recursively all polynomials of degree < m. 

Examples. The case q = 2. We get a sequence of polynomials that we will 
call the central Bernoulli polynomials (of the first kind) : 



The generating function reads 

oo r>II 



2^ fc! 



fc! \sinht 

fe=0 



t 



e xt . (6.8) 



For x = one obtains the numbers B^ . From the generating function it 
emerges that Bl 1 = for fc odd. The first central Bernoulli numbers B^l are: 

1 _I JL _31 127 2555 t 
A ' 3 ' 15 ' 21 ' 15 ' 33 ' eLL- 



The polynomials (|6.7|l can be directly expressed in terms of the classical 
Bernoulli polynomials: 



B» (2x - 1) 



Other properties are now briefly discussed. From the polynomial expansion 
Theorem we get 



*»f. (x) = — u (6-9) 



fc=0 

For a = 1, since 



e* - e"* 



fc-i 



v = n (fc - 1)!5 7/ (n - 1, k - 1) , (6.10) 



we have: 



Kla (a) = <a (0) + £ (n - 1, * - 1) . (6.11) 



fc 

k=l 
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For q = 3, the corresponding delta operator is given by formula (|4. lUp> . with the 
generating function 

°° fill I ( \ f xt 

2^ fc! e 2 t _ 2e t + e-*- 1 j 

fe=0 

For the other cases, we get similar results. For instance, 

f = ^ (6 13) 



k\ 

k=0 

°° f>VII f \ _ fp xt 

V [ -^-t k = - (6 14) 

k=0 

and so on. As an example, the first polynomials i?„ (x) are given by 
Bi u (x) = 1, B[" (x) = x + | B{ u (x) = x 2 + 3x+^, 

B\ u (x) = x 3 + % -x 2 + y X + 39 ' 
B{ u (x) = x 4 + 6x 3 + 37x 2 + I56x + HZH 

Bl" (x) =x 5 + -x 4 + —x 3 + 390x 2 + —x + 3365, 
5 w 2 3 6 

B\" (x) = x 6 + 9x 5 + ±fV + 780x 3 + ^-x 2 + 20190x + U6 ™ 7 
D v ' 2 2 42 

Remark 6. The polynomial sequences generated via the formulae H6.12J) - I|6.14JI 

and many others obtained according to the previous construction possess the 
remarkable property that c p _i = 1 mod p. Therefore, they satisfy the hypothe- 
ses of Theorem ^ and provide examples of polynomial sequences verifying the 
Almkvist— Meurman congruence. 

The relation (|6.11|) between higher order Stirling-type numbers, Bernoulli- 
type numbers and Bernoulli-type polynomials can be generalized as follows: 

n 

B? ha (x) = Bl a (0) + T S " (» - 1. k ~ 1) • ( 6 - 15 ) 
fc=i 

The recurrence relation for the polynomials B^ a (x) is formally derived from 
the defining relation (|6.1I) : 

f>e^. Q (x) = ( J ) tx n = nBl_ ha _ x (x) . (6.16) 
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Therefore, the difference equation solved by the Bernoulli-type polynomials is, 
for every q 

m 

E a l B n,a ( X + 3) = nB l-l,a-l (*) , (6-17) 
3=1 

where the satisfy the constraints (|4.2|) . 1)4. (|4.4|l . When q = 2, we obtain 
(aj + 1) - (x - 1) - nBiUa-i (*) = 0. (6.18) 

The recurrence relation (|2.5|) valid in the case of universal Bernoulli polynomials 
reads: 

Bl +ha (x)=(x-^jBl a (x), (6.19) 

where 

Using the relations [t,x] — 1 and B% a — g (i) -1 1™, which is a consequence of 
eq. p.lOfl . we get 

(n + 1) S* (z) = f xi + 1 - o ^'T^^t" 1 ^ S ^ (*) ■ ( 6 ' 2 °) 

From the previous equation we derive the formula expressing the Bernoulli-type 
polynomials of order a + 1 in terms of those of order a: 



E Mia+i (Z + k) = (l - ^) Bl >a (x) + ™B£_ li0 (x) . (6.21) 



Observe that, in accordance with eq. (|4.15() . the polynomials 

Pn (x)= x ( - f x n -' - xB*_ lt na (x) (6.22) 

represent the basic sequence associated with the operator Q = (4r) Q D. Con- 
sequently, for a = l, we deduce the relation expressing the basic sequence x^ in 
terms of Bernoulli-type polynomials: 

(x n y = xBl_ hn (x) . (6.23) 

In perfect analogy with the case of the classical Bernoulli polynomials, an 
interesting connection between the Bernoulli numbers and the Stirling numbers 
of the same order holds. Indeed: 

s" (n,r) = ±{t r \(x n y) = ±(t°\f(x n y) 



^ Bn- r ,n+l (-0 : 



(6.24) 
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or, equivalently, 



Therefore, 



Analogously, 



* 9 (n,r) = ( U r _ \ )i?«_ r , n (0). (6.25) 



1 / (T,k=i a ke 



kt\ r 



S* („, r) = £- ( ( ^feLT^ ) | fV l ) = ( " ) B^U-r (0) , (6.26) 



E " ) fl «-r,-r(0)(^) 9 - (6-27) 



and 



7 Sequences of integer numbers 

In this section, sequences of integers are constructed by exploiting the properties 
of the universal Bernoulli numbers previously discussed. 

Lemma 5. Consider a sequence of the form 

1 jr§t\ (7.D 



Gi (t) G 2 (t) ^ k\ 

where G\ (t) and G 2 (t) are formal group exponentials, defined as in formula 
fj.6j) . Assume that c n 6 N and c p _i = 1 modp, when p > 2 is a prime number. 
Then {Nk}keN is a sequence of integers. 

— l 

Proof. It suffices to observe that the Bernoulli-type numbers Bk and 

Bk associated with the formal group exponentials G\ (t) and G 2 it) , under the 
previous assumptions must satisfy the classical Claussen-von Staudt congruence 
for k even, as a consequence of Clarke's universal congruence (|2.7H . For k odd 

these numbers are also integers, due to 1)2. 8[l . It follows that the difference Bk — 

Bk for any k is an integer. 

Corollary 2. Assume that the hypotheses of Lemma\^hold. Then any sequence 
of polynomials of the form 



t t 



Gi (i) G 2 (t) 



fc=0 

is an Appell sequence with integer coefficients. 



e *t = £ 3M t k (7.2) 
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For instance, the sequences associated with l|6.13[) and (|6.14|) satisfy the 
Claussen - von Staudt congruence in the strong sense and can be used to con- 
struct integer sequences of numbers and polynomials. Here are reported some 
representative of generating functions of integer sequences which can be ob- 
tained from the previous considerations. 

a) 

t(l + e«) 
1 + e* - e 2t ' 

The sequence generated is 

6, 114, 2901, 95436, 3894815,... 

b) 



t (1 + 2 cosh t + 4 cosh 2t - 6 sinh t) 



( - 3 + 4 cosh t) (2 + cosh t - cosh 2t - sinh t + sinh 2t + 2 sinh 3t ) ' 
The sequence generated is 

-7, 61,-642, 10127,-207110, 5001663,... 

c) 

te"f sec/i(f ) [4 + e*(l + 2e*(-l + e*))] 
2(-3 + 4 cosh t) (1 + (-2 + 4 cosh t) sinh t)) ' 
The sequence generated is 



-5, 29,-150, 1279,-17770, 268647,. 



8 New Bernoulli— type polynomials of the sec- 
ond kind 

By analogy with the case of classical Bernoulli polynomials of the second kind 
( |39p. we introduce the following sequences of polynomials. 

Definition 14. The higher-order Bernoulli-type polynomials of the second kind 
are the polynomials defined by 

bl(x) = J q (t)(x)l, neN, ?£N, (8.1) 

where J q (t) is the operator \b\ fy) . 

These polynomials represent a Sheffer sequence for the operator A 9 of eq. 
6U, (IP) : 

A q bl(x)=nbl_ 1 (x), (8.2) 
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and therefore they satisfy the identity 



bl(x + y) = J2[k K («)«-*. (8-3) 



k=0 

which relate them with the higher factorial polynomials. In particular, for y = 0, 
we get 

%(«)=E( I )««_ fe - (8-4) 

k=Q ^ ' 

From eq. I|8.2|l we deduce the difference equation satisfied by polynomials H8.1(l : 

J2a k b q n (x + k)=nb q n _ 1 (x). (8.5) 

k 

The case g = 1 just reproduces the standard Bernoulli polynomials of the second 

kind. The generating function can be obtained in some specific cases. Let us 

j j t _ —t 

briefly discuss the second order case. Since the operator J = 2— =f — acts as 
follows: 



f 



x + l 



J P( x ) = 2 / P ( u ) du > 



x-1 



we deduce an explicit expression for the central Bernoulli polynomials of the 
second kind: 



x+l 



b I n I (x) = J II x!! = - ^ i (u) n du. (8.6) 
The generating function is 

> = — 7 . (f + yi + t 2 ) . 8.7 

^ fc! log (t + vTTF) V J 

For x = we get 

^ & J fc J (Q)* fc _ t 

*l io g (t + VTTF)- l ' j 

Let us derive the recurrence relation satisfied by the polynomials 6„ (x) . From 
the defining eq. (|8.6I) we get 

K 1 (*) = (^^) (*)» = (8-9) 

Integrating we obtain 

# (x) - # (0) = n I* (utiUdu. (8.10) 
Jo 

The difference equation (|8.5(l reduces to 

b'J (x + l) = bi 1 (x-l) + 2n&£ 1 (») • (8.11) 
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Let us compute the numbers 6„ (0) and 6„ (1). For x = 1, we get 

# (2) = # (0) + Inb^U (1) - bi 1 (0) + 2n [ (uY^du, (8.12) 



so that 

Jo 

and ^ 

Jo Jo 

Notice the analogy with the classical Bernoulli numbers of the second kind, 
defined by gHj 

bn (0) = (^—^ | (x) n J> = J 1 (u) n du. (8.14) 

There is a connection between Bernoulli-type polynomials of second kind and 
generalized Stirling numbers of the first kind. Since 



n 1 

we get 



fc! 

fc=0 



K (*) = bl (0) +£V( n -l,k- 1) x fe . (8.15) 



/fc 
fe=i 



Many other properties and identities, that we will not discuss here, can be 
derived using operator techniques. 



9 New Euler— type polynomials and numbers 

In this Section, a new class of polynomial sequence of Appell type is discussed. 
For many aspects, it can be considered to be a natural generalization of the Euler 
polynomials. The Definition proposed in this paper is different with respect to 
that in 

We recall that the Euler polynomials of order a, which will be denoted by 
the symbol (x) , are the Appell sequence associated with the operator 




a ^ 0. 



Consequently 

K (x) = (l + x n . (9.1) 

Therefore, we can propose the following 



23 



Definition 15. The Euler-type polynomial sequences are defined by the relation 

(x) = (l + ^j x n , (9.2) 

where a,u> £ K. 

Of course, this Definition is equivalent to that proposed in the introduction 
in terms of a generating function: 

y H W t k =1 + ^1) e\ (9.3) 

For the sake of clarity, we will omit the superscript uj. The Euler-type 
polynomials satisfy the Appell property 

DET (x) = nE%\ (x) 

and the Appell binomial identity 



i — n V / 



x n ~\ (9.4) 

fc=o ^ ' 

As in the standard case, this enables to define Euler-type polynomials in terms 
of Eulcr numbers: 

£r(z) = E( I W a (°)^ fc - (9-5) 

k=0 ^ ' 

It also emerges that for every q the polynomials H9.3|) are an Appell cross- 
sequence: 

Ek a+b ^ + v) = jz{l ) E u a (y) (*) ■ (9-6) 

fc=0 ^ ' 
We also have a recurrence formula : 

K'ti (*) = ( x ~ ^||) E T (*) . (9-7) 

where g q (t) = {l J- 
There exists an interesting analogue of the famous Boole summation formula. 
Since, for any function h (t) 

Mt) _£<*wi*p» (hMV, 

fe=o ' V w / 

then, choosing /i (t) = exp (yi) we get 

^M = £^^(l + (9-8) 
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Once applied to a polynomial p (x) , in the case a — 1 we obtain for any q the 
expansion 

^+ y ) = E^^(i + ^) Q p« (°)- 

fc=0 ' ^ u ' 
Let us first consider the operator (for q = 2) 

g II (t)= (l + l^flj. (9.9) 

Correspondingly, 



' 1 + < ~J~) (9 ' 10) 



First we observe that, since 



1 \ sr^f-a\ I e — e x 



the polynomials l|9.10|l are of the form 

«" « - £ (7) ( f ^ fZ! )^" = t £ (7) iS s " <*■*> 

j=0 v J 7 v 7 fc=0j"=0 v J 7 

(9.12) 

From this equation and from the Definition (|9.1Q(I we derive the difference equa- 
tion 

(x + 1) - i^' a (i - 1) = ui^- 1 (a:) . (9.13) 
The analogue of the Newton expansion is given by the formula 

n n , s 

which is easily derived from eq. (|9.12|) and from the formula 



f —t \ 
e - e 



x"=^^(n,fe)(i) fe x^ fc . (9.15) 
Finally, from the recurrence formula l|9.7|) (for q = 2) we have: 

providing another difference equation satisfied by E^ I,a (x): 

^n+l (s) = ^^ 7 /' Q (*) - - (^' a+1 (*+!) + ^' a+1 (x 1)) . (9.16) 
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